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We observe strong, non-trivial cross-correlations between 
the edge states found in the interior and the exterior of mag- 
netic quantum billiards. Our analysis is based on a novel 
definition of the edge state spectral density which is rigorous, 
practical and semiclassically accessible. 



One of the main goals in the field of "quantum chaos" 
is to link the auto-correlations found in a quantum spec- 
trum to the periodic orbits of the classical problem [Q. 
Here, we extend this study and investigate whether cross- 
correlations exist between quantum systems which are 
different but related by their classical dynamics. We de- 
velop this idea for magnetic quantum billiards which 
often serve to model semiconductor quantum dots |^]. 

Magnetic billiards consist of a charged particle moving 
ballistically in a compact domain in the plane subject 
to a homogeneous magnetic field. The quantum wave 
function is required to vanish at the billiard boundary 
while the impinging classical particle is reflected specu- 
larly 1^ ||. The boundary defines also a complementary 
problem - an anti-dot - where the particle is confined 
to the exterior, and is scattered at the billiard boundary. 
Although the exterior domain is unbounded, its spectrum 
is discrete Is it possible to relate the energy levels of 
the dot to those of the anti-dot? We shall show that 
there exists an intimate, non-trivial connection between 
the spectra of the interior and the exterior problem. It 
is the quantum manifestation of a duality in the classical 
dynamics. 

The classical interior- exterior duality is illustrated in 
Fig. |l|(a): Since a (periodic) orbit consists of arcs of con- 
stant curvature, one can construct a dual orbit in the 
complementary domain by completing the arcs to cir- 
cles. Any skipping trajectory meets with a dual one un- 
der rather general conditions - if every circle of cyclotron 
radius p intersects the boundary at most twice. Pairs of 
dual periodic orbits have the same stability and their 
actions add up to an integer multiple of the action of a 
cyclotron orbit. On semiclassical grounds one may there- 
fore expect the correlation between the interior and the 
exterior motion to carry over to the quantum spectrum. 
This is also corroborated by the existence of pairs of in- 
terior and exterior quantum eigenstates which match up 
well, cf Fig. |l|(b), although their energies differ. 

The semiclassical analysis is complicated by the fact 
that in the exterior each Landau level is an accumula- 



tion point for an infinite series of energies. The respective 
states - the bulk states - correspond to unperturbed cy- 
clotron motion. Also in the interior one may find (a finite 
number of) bulk states if p permits complete cyclotron or- 
bits to fit into the domain. The eigenfunctions which cor- 
respond to the skipping trajectories, on the other hand, 
are called edge states. Clearly, a possible correlation is to 
be expected only between these non-trivial exterior and 
interior states. 
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FIG. 1. (a) A pair of dual periodic trajectories and (b) a 
pair of correlated eigenfunctions found in the interior and the 
exterior ellipse billiard (superimposed). The billiard bound- 
ary is indicated by a dashed line. 

Although the notion of edge states is intuitively clear 
and often used (e.g. in the context of the quantum Hall 
effect jl^), we are not aware of a general quantitative 
definition in the literature (however see |11|). Therefore, 
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the purpose of this letter is twofold. First, we propose 
a definition for the spectral density of edge states which 
provides a meaningful characterization and applies in the 
quantum and in the semiclassical regime. Only with this 
can we then establish the existence of a pairwise relation 
between the edge states of the interior and the exterior. 

A definition for edge states should take into account 
that a clear separation into edge and bulk occurs only 
in the semiclassical limit 6 — > 0. Here, we express 
the quantum scale in terms of the magnetic length b = 
y^2h/{mujc) (with m the mass and ujc the cyclotron fre- 
quency). At finite values of b the states may be of an 
intermediate type. We propose to quantify the transi- 
tion by attributing a weight > to each eigenstate 
ipi (of energy Vi) which gives a measure of the degree to 
which ipi has the character of an edge state. The spectral 
density of edge states in either the interior or the exterior 
is then defined as 
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Here, we scale the energy E by the spacing between Lan- 
dau levels, V = E/{hwc) ~ 1^- Any reasonable def- 
inition of the weights Wi must suppress the bulk states 
by exponentially small values, such that the mean edge 
density dedge is well-defined in the exterior and equal to 
the interior one, to leading order. In the semiclassical 
limit it should match our notion of edge states admitting 
a trace formula which involves only the skipping trajec- 
tories. Moreover, we shall demand dcdgo to coincide with 
the unweighted interior mean density if the cyclotron ra- 
dius is large enough to prevent bulk states. 

To motivate our definition of the weights Wi, consider 
the scaled magnetization A4 of the interior billiard, a sum 
over the scaled magnetic moments 
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The scaled magnetization density m(i/) can be expressed 
by the derivatives of the spectral counting function 
N(z/; b) = Q{v — Vi{b)) with respect to b^ and 
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This is verified by replacing the energies Vi in (||) by 
the expectation values of the Hamiltonian. The scaled 
magnetization density may be obtained from the conven- 
tional one by multiplication with the field strength B. It 
exhibits a natural partitioning into a bulk part and an 
edge part since it complies with the scaling properties of 
the system: The scaled magnetic moment of a Landau 
state is —v. Hence, the second part of (0), 
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9N 
dv 



]{-Vi)5{v - Vi) , 



(4) 



attributes the full diamagnetic response of a Landau state 
to each state V'i- We call it the bulk magnetization den- 
sity. It follows that the remaining part, the edge magne- 
tization density 
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assigns the positive excess magnetic moments induced by 
the presence of a billiard boundary. Its mean value 
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follows from the mean number of states in a magnetic 
billiard with area £^ and circumference 1l2l , 
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Note that lUcdgo cancels the mean bulk magnetization 
density exactly: There is no orbital magnetism apart 
from the quantum fluctuations. Hence, mcdgc character- 
izes those few (edge) states which carry a finite current 
along the boundary, balancing the bulk magnetization 
due to their large positive magnetic moments. 

The edge magnetization is well-defined in the exterior 
as well. There, it is negative with the mean like (^ but 
for a minus sign in front of the area term. This suggests to 
define the edge state density as dodge (i') = ±modgc(i^)/i^, 
with the lower sign for the exterior problem. The corre- 
sponding weights 
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are easily obtained as the derivative of the eigenenergies 
taken at fixed p. This definition satisfies the conditions 
formulated above. In particular, the weights are expo- 
nentially small for bulk states since the Landau energies 
V = N + \, N ^^Q, are independent of 6 

The semiclassical edge state density is derived by in- 
serting the trace formula for N(i/) ||l4|,|l^ into (||). One 
obtains a sum over all skipping periodic orbits for the 
fluctuating part dedge = dedge — dedge- It differs from 
the semiclassical expression of the unweighted spectral 
density only by a factor 



(9) 



attributed individually to each periodic orbit contribu- 
tion. This classical weight is determined by the area ^ 
enclosed by the trajectory 7 and by its length Jf^. The 
weights approach zero as the skipping orbits are further 
detached from the boundary. Hence, the classical weights 
(^ smoothly suppress the bulk contributions to the semi- 
classical spectral density. 

Fig. ^ shows (a) the quantum weights Wi against their 
energies Vi [p| and (b) the phase space distribution of 
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the classical weights (for the interior ellipse billiard with 
area = -k and eccentricity 0.8 at 6 = 0.1; the shade 
in (b) gives the probability measure for finding a tra- 
jectory with weight w-^.) In Fig. |2|(a) one observes how 
the Wi distinguish the relevant edge states from the bulk. 
Characterized by vanishingly small Wi, the bulk states 
accumulate at the Landau levels, v = + i, with se- 
quences of transitional states connecting to the large edge 
weights. The latter distribute in structures reproduced 
by the classical weig hts (|) (which are due to bifurcat- 
ing regular islands in phase space). In the exterior, the 
segregation into edge and bulk is even more distinct (not 
shown). 
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FIG. 2. (a) Weighted edge state spectrum compared to (b) 
the phase space distribution of the classical weights for the 
interior ellipse. The quantum weights (|^ segregate the edge 
states from the bulk and mimic the structures in the distri- 
bution of classical weights (|^. 

In order to unravel the relation between the interior 
and the exterior spectra we consider the cross-correlator 
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defined for fixed reference energy and magnetic length 
h{j. The functions g and h are normalised Gaussians. 
g regularises the integral over the pair distribution 



^cdgc('^! '^odgc ^ restricts the 6^ integration 

to a range where a linear expansion of Vn{b^) near 6q is 
allowed. Then, C(z/o) can be written as a double sum 
over the interior and the exterior edge spectra obtained 
at one fixed value bo of the magnetic length: 
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The primes label the exterior energies and weights. Since 
the width of g is taken small, ag ^ (^cdgo)""'^, only those 
pairs of interior and exterior energies contribute whose 
distances to vq scaled by the the respective quantum 
weights are approximately equal. The prefactor in (|TT 
ensures that only pairs of edge energies contribute. 

When evaluating C(i'o) semiclassically, the integration 
over in (|l0|) selects those pairs in the sum over inte- 
rior and exterior orbits 7 and 7' which satisfy w-yJi'-y ~ 
Wy'J^^i, a relation fulfilled by the dual pairs. Restricting 
the summation to the latter is tantamount to the "diag- 
onal approximation" p8| . The actions complement each 
other to 27ri'on-y where 2'kvq is the scaled action of a cy- 
clotron orbit, and is the number of reflections. We 
obtain 
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with g the Fourier transform of g. Here, rimin is the 
minimal number of reflections needed for a periodic orbit 
at given p and 
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a sum over the set r„ of dual orbits with n reflections. 
It involves the classical weights w~f (^) and the stability 
amplitudes a-y |14| of the unweighted spectral density. 

From its semiclassical form (|l2|) the correlator is ex- 
pected to be appreciably different from zero only at ener- 
gies where the cosine terms are stationary. Hence, C(fo) 
must exhibit peaks &i vq — N + ^. Its Fourier transform, 
on the other hand, should be peaked at the integer values 
starting from nmin- 

A numerical verification of these predictions is pre- 
sented in Fig. ^ for the spectra of the ellipse billiard at 
&o = 0.1 ^. In Fig. ^(a) one observes that the cross 
correlation function is strongly fluctuating but displays 
pronounced spikes at the expected energies. They are a 
clear signature of a correlation between the interior and 
exterior edge states. The peaks in the Fourier transform. 
Fig. H(b) , are positioned at integer values of t which start 
at = 4, as expected for the (desymmetrized) ellipse. 
They expose clearly the classical duality as the origin of 
the cross-correlations. 
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As a test we restricted the sum ( 11 ) to pairs taken from 
the two different symmetry classes of the eUipse. This 
erases the peaks in C{i'q) as one expects semiclassically. 
On the other hand, removing the bulk states by imposing 
a threshold on the Wi does not change Fig. 



state spacings, Vi = 31.42554, v'^ = 31.61696, the differ- 
ence between the two sides of (14) with Wi = 17.23/z/i, 




FIG. 3. (a) Cross-correlation function ( |ll| ) for the ellip- 
tic billiard. The pronounced spikes indicate the existence 
of pairwise correlations between interior and exterior edge 
states, (b) Fourier transform of (|l^). The peaks at integers 
starting from 4 prove the classical origin of the correlations 
(0-9 = 0.001). 

The spikes of C(i'o) imply the existence of a pairwise 
relation between the interior and exterior edge states: For 
each correlated pair of interior and exterior edge energies, 
i^i and i/'j , there exists a Landau level I'o = + 5 such 
that the distances - scaled individually by the quantum 
weights Wi and w'j - are approximately equal. 
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This follows immediately from (|11|), where we took the 
width of g to be small on the quantum scale. The relation 
( |l^ ) shows that the interior and exterior edge spectra are 
intimately connected in the semiclassical limit. Note in 
particular the vital role played by the quantum weights 
without which the correlations would not be observ- 
able. 

Equation ( p^ allows to spot single pairs of correlated 
states in the spectrum, and Fig. ^b) gives an outstand- 
ing illustration (with the shade proportional to IV'P): 
The wave functions are clearly localized along the sta- 
ble dual periodic orbits drawn in Fig. |l](a). Although 
the respective energies are separated by 20 mean edge 
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— 26.76/j^j, and = 31 is approximately one tenth 
of the mean level spacing scaled by the mean weight. The 
correlation of pairs of chaotic wave functions cannot be 
verified as easily by visual inspection, but they exhibit a 
large overlap of their normal derivatives at the boundary. 

Let us finally emphasize that the correlations between 
edge states of the interior and the exterior do not per- 
mit us to derive one spectrum from the other - even 
in the semiclassical limit - since the respective Landau 
level N and the complementary quantum weights are not 
known a priori. Nonetheless, they quantify a deep in- 
terconnection of the spectra which is generated by the 
classical duality, as initially conjectured from Fig. ^. 
Similar cross-correlations may be expected between non- 
magnetic systems as well, e.g. in complementary billiards 
on the sphere. 
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